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Abstract 

A phenomenological quasiparticle model, featuring dynamically generated self-energies of exci- 
tation modes, successfully describes lattice QCD results relevant for the QCD equation of state 
and related quantities both at zero and non-zero net baryon density. Here, this model is extended 
to study bulk and shear viscosities of the gluon-plasma within an effective kinetic theory ap- 
proach. In this way, the compatibility of the employed quasiparticle ansatz with the apparent low 
viscosities of the strongly coupled deconfined gluonic medium is shown. 



1. Introduction 

Within the past years, an immense effort has been put into revealing the nature and proper- 
ties of deconfined strongly interacting matter both theoretically and experimentally in relativistic 
heavy-ion collisions. The aim is to understand characteristics of QCD matter such as its equation 
of state (EoS), collective behaviour or transport properties which have a wide range of implica- 
tions in cosmology and astrophysics. While the equation of state describes the system in thermal 
equilibrium, the transport coefficients, such as viscosities, characterize the system's ability to 
relax from nonequilibrium towards equilibrium. 

The success of perfect fluid hydrodynamic calculations in describing the collective flow ob- 
served in heavy-ion collisions yj, |2, l3D suggests small viscosities of the created medium which 
has also lead to viewing it as being strongly coupled |4]. Nonetheless, in particular the shear vis- 
cosity cannot be arbitrarily small due to unitarity Jj|]. Besides, in weak coupling certain relations 
hold among different transport coefficients, cf. e. g. 0]. 

Apart from the first-principle calculations of bulk (0 and shear (rj) viscosities by means 
of lattice QCD Q H a variety of other approaches was proposed. Among these the rigorous 
perturbative calculations starting either from the Boltzmann equation ]9l LL0|] or from the Kubo 
formalism full have to be mentioned. Besides, considerations employing different spectral func- 
tions lfl2ll or numerical transport calculations of the Boltzmann equation El have been made. 

Here, we adress viscosities by viewing the gluon-plasma as composed of quasiparticle ex- 
citations. The underlying quasiparticle model (QPM) was successfully tested to describe lattice 



QCD results of the EoS 11141 11511 . In the QPM, gluon and quark quasiparticles obey dispersion 
relations, where the entering self-energies n, in general, depend on temperature T and chem- 
ical potential yt both explicitly and also implicitly via a phenomenological effective coupling 
We extend this picture to nonequilibrium systems by means of an effective kinetic 
theory for gluon quasiparticles, i. e. // = in the following. 
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2. Effective kinetic theory for quasiparticle excitations 



For a system in nonequilibrium, the gluon quasiparticle dispersion relation no longer depends 
on an uniquely defined temperature but becomes space-time dependent, E(x) = -\Jp 2 + H(x). 
In this case, the space-time dependence of the distribution function b(x, p) is governed by the 
Boltzmann equation 

lp a (x)d a + yfn(x)F a (x) jJL^ J b(x, p) = Db{x, p) = C[b(x, p)], ( 1 ) 

where C[b(x, p)] is the collision term. The force F a = d a Tl/ (2 VfT) satisfies p a F a = such 
that the spatial gradient of the self-energy acts as an external force changing the momenta of 
the quasiparticles between collisions. The scattering interaction conserves locally energy and 
momentum which results in a vanishing collision term when multiplied by p v and integrated 
over three-momentum p. Hence, an energy-momentum tensor 

T " v ( x ) = d f - d l P - pf(x)p v (x)b(x,p) + g^B(Tl(x)) (2) 
J (27TyE(x) 

can be defined which, as a consequence of the Boltzmann equation, obeys energy-momentum 
conservation dfjT^ix) = under the condition 

dB 1 f d 3 p 

TJPfTT = -~q(x) , q(x) = d b(x,p) , (3) 

oYl(x) 2 J {27iyE(x) 

for the mean field B. Here, q(x) is an auxiliary field = diagi 1,-1,-1,-1) and d is the 

number of degrees of freedom. As a result of Eq. (0), the space-time dependence of the self- 
energy is determined by the auxiliary field. Thus, the Liouville operator T> on the left hand side 
of Eq. (Q} is a functional of the distribution function b(x, p) itself. 

In thermal equilibrium, characterized by a local distribution function b°(x,p) = (e p °""l T - 
where the fluid four- velocity u a satisfies u a u a = 1, one has to demand that Yl{q(x))\p = 
n(T) to recover equilibrium results from the energy-momentum tensor ©. In fact, by comparing 
Eq. (01 evaluated for b°(x, p) in the local rest frame with T*Z. = eu f 'u v -P(g f ' y -u f 'u y ), the quasipar- 
ticle model expressions for energy density e and pressure P are recovered, cf. 114, 11611 . Moreover, 
Eq. (O represents the nonequilibrium generalization of the stationarity condition ensuring ther- 
modynamic self-consistency and also implying that the principles of statistical mechanics assure 
the physical meaning of b°(x,p) (lib . In this way, T^ix) in Eq. (|2]l represents the general form 
for an isotropic fluid composed of quasiparticle excitations including only g^ v and u 1 " which, in 
addition, satisfies thermodynamic self-consistency in equilibrium. 

3. Bulk and shear viscosities 

The calculation of transport coefficients from T^ix) as the first-order corrections to thermal 
equilibrium is straightforward. Assuming small deviations from equilibrium, b(x, p) = b°(x, p) + 
6b(x,p) with 6b <K b°, T^ v {x) can be decomposed into = T^[b°] + 5T^ v [Sb]. Expanding 
B(Tl(q(x))) in terms of small deviations from its equilibrium value and approximating E(x) to 
lowest order by E — y/p 2 + n.(T), ST^ to lowest order in 5b can be written as 



Assuming that collisions always result in an exponentially fast restoration of local equilib- 
rium with the relaxation time r, all the complexity of C[b] is encoded in r Correspondingly, 
6b can be approximated by 6b = —C[b]/T, which in turn can be expressed by the Boltzmann 
equation ([]]). Thus 6b in Eq. (01) is a functional of D acting on b reading to lowest order 

6b(x,p) = -I [p a d a - UmT)^b\x,p) . (5) 

Considering gluonic quasiparticles with Tl(T) = jT 2 G 2 (T), we note first that for recov- 
ering the equilibrium QPM one needs to identify 11(g) = qlifio \og[M -\/q - T s }/T c ] 2 ), where 
q\ b o = (Nq)\y> = r 2 and/3 = 11/(8tt 2 ). Furthermore, one finds q (dll/dq)\ b0 = T 2 (dYl(T)ldT). 
When evaluating Eq. © supplemented by Eq. ©, it is appropriate to replace the convective 
derivatives of T and it by spatial gradients. In line with the Chapman-Enskog strategy, the space- 
time dependence of b is assumed to be determined only by T and u^, constituting b°, as well as 
their gradients. Then, the conservation equations are given in first approximation by the equa- 
tions of motion and energy of an ideal fluid. To make the decomposition of b(x, p) unique, we 
employ the Landau-Lifshitz condition lfl9ll . i. e. UyT 1 ™ = u v T?Z = eii 12 , implying ST 00 — in the 
local rest frame. 

For the considered gluonic system, the bulk and shear viscosities are the only independent 
transport coefficients characterizing the fluid. In case of small deviations from equilibrium, 
they are obtained from the spatial part of the nonequilibrium energy-momentum tensor ST'i = 
-£6 l idkii k + TjW'-i as coefficients of the scalar and the traceless part W-* = d l u' + d'u 1 - 2 :g'-'di l u k . 
In the local rest frame one finds 



dE 

E - T — 

dT 



From Eqs. (O and Qi it is clear, that only £ is significantly influenced by the medium-dependent 



d f d : 
15TJ (2^r 



quasiparticle dispersion relations E, cf. also 12011 . while in rj the mean field contributions vanish. 



To quantify r\, for instance, the parameters of the effective coupling, T s and A/T c , are first 
adjusted to lattice QCD results of the EoS, for example, of the scaled interaction measure (e - 
3P)/T 4 (left panel of Fig. [I}- In addition, for r required in Eq. (0, we employ an ansatz inspired 
by previous work [22], t~' = a rj /(32n 2 )TG 4 log(a ;; 7r/G 2 ), where the QCD running coupling is 
replaced by our effective coupling G 2 (T). The corresponding result for the shear viscosity to 
entropy density ratio rj/s is depicted in Fig.Q](right panel). 



4. Conclusion 

We discuss the bulk and shear viscosities by means of an effective kinetic theory for gluon 
quasiparticle excitations in the relaxation time approximation. Thermodynamic self-consistency 
necessary in such a quasiparticle description turns out to be solely a consequence of energy- 
momentum conservation. Our numerical results for rj/s are in agreement with available lattice 
QCD results J3>Sl and with other approaches Moreover, rj/s exhibits the expected mini- 



mum close to T c similar to classical fluids 02311 . which here is mostly driven by r. Contrary to 



the popular view, where a large quasiparticle mean free path implies large 77/s 024I1 . our results 
suggest that a quasiparticle description is still admitted for the strongly coupled gluon-plasma. 
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Figure 1: Left: Comparison of QPM results for the scaled interaction measure with lattice QCD results for pure 
SU(3) (2lll (boxes). The adjusted model parameters read d = 16, T s = 0.53 T c , A = 2.88, where in line with lattice 
QCD results we set T c = 271 MeV. Right: Corresponding QPM result for T]/s, employing a, ; = 6.8 in the ansatz for the 
relaxation time r, compared with lattice QCD results (boxes from [7], diamonds and triangles from [8]). In addition, the 
unitarity limit [5] rj/s = 1/(4jt) (dotted curve) is depicted. 
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